Characterization of addressability by simultaneous randomized benchmarking 
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The control and handling of errors arising from cross-talk and unwanted interactions in multi- 
qubit systems is an important issue in quantum information processing architectures. We introduce 
a benchmarking protocol that provides information about the amount of addressabihty present in the 
system and implement it on coupled superconducting qubits. The protocol consists of randomized 
benchmarking experiments run both individually and simultaneously on pairs of qubits. A relevant 
figure of merit for the addressability is then related to the differences in the measured average gate 
fidelities in the two experiments. We present results from two similar samples with differing cross- 
talk and unwanted qubit-qubit interactions. The results agree with predictions based on simple 
models of the classical cross-talk and Stark shifts. 

PACS numbers: 03.67.Ac, 42.50.Pq, 85.25.-j 



The ability to selectively address one qubit or subsys- 
tem from among many in a quantum register is a key 
prerequisite for a scalable quantum computing architec- 
ture. This addressability can be lost either due to con- 
trol fields that target one subsystem influencing the other 
neighboring subsystems (classical cross-talk errors) , or by 
unwanted quantum interactions between the target sub- 
system and the other subsystems. Simple estimates for 
the classical cross-talk errors can be obtained by mea- 
suring the residual Rabi rate on un-targeted qubits [ij. 
Extensions to capture all the relevant error rates are nec- 
essary. This issue is particularly important because tra- 
ditional models for fault-tolerant quantum computation 
assume uncorrelated errors and therefore these address- 
ability errors need to be characterized and subsequently 
minimized. 

The experimental demonstration of errors of the order 
required for fault tolerant quantum computation 
presents a formidable challenge, as it is typically limited 
by errors in state preparation and measurement (SPAM) . 
Indeed, quantum process tomography (QPT) one 
of the standard methods for characterizing a quantum 
process, is particularly sensitive to SPAM errors. Addi- 
tionally, complete quantum-gate characterization rapidly 
becomes experimentally intractable due to the exponen- 
tially large Hilbert space. 

An alternative is randomized benchmarking (RB) 0- 
[loj . which has been used to characterize single qubit 
gates in liquid-state NMR trapped atomic ions 

and superconducting qubits [l5l - [l7| . RB is 



specifically tailored to compensate for SPAM errors by 
considering only the exponential decay of sequences of 
random gates. This comes at the cost of only obtaining 
information about the average gate error over the Clif- 
ford group, although some alternative approaches have 
been recently devised for extending RB to estimate the 
error of a single, particular gate 18. T^. 



In this Letter, we present and experimentally imple- 
ment a method for characterizing the amount of address- 
ability between two subsystems. We typically specialize 
to qubits, but the method is general and, for example, 
handles subsystems with logical qubits. First, we per- 
form RB on each of the subsystems independently, while 
leaving the other unperturbed. Then we perform (sub- 
system local) RB on both subsystems simultaneously, 
and compare the error rates. The comparison of these 
two can indicate a parametric dependence of the error 
rate of a subsystem on the operations in another subsys- 
tem. Furthermore, we show how the twirling protocol 
proposed in Ref. [20[ can be modified with this RB pro- 
tocol to detect some forms of spatial correlations. These 
methods utilize many of the techniques of Ref. 0] and 
are similarly scalable and less dependent on SPAM errors 
than QPT. 

We start by defining a general randomized benchmark- 
ing protocol. We choose a sequence of m-|-l unitary gates 
where the first m gates are chosen uniformly at random 
from a group Q = {lAi\ and the (m-fl)"^ gate is chosen to 
be the inverse of the composition of the first m random 
gates. Provided Q is at most the Clifford group, this undo 
gate can be found efficiently by the Gottesman-Knill the- 
orem [i^l- Assuming each gate Ui has some associated 
error, represented by A,;(/9), the sequence of gates is mod- 
eled by 



5i„ 



A,- 



(1) 



where o represents compositions, im is the m-tuple 
(ii, im) describing the sequence, and the recovery gate 
^im+i is uniquely determined by im- Since {Ui} forms a 
group, the sequence can be rewritten as 



= A, 



0-1 



o A, , oU::, 



(2) 



where Zi;. is another clement of Q. 
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For each sequence the overlap 'Ti[ESi^{p)] between 
-E, an operator representing the measurement (with er- 
rors), and Si^{p), the final state, is measured. Av- 
eraging this overlap over K independent sequences of 
length m gives an estimate of the average sequence fi- 
delity Fscq{m, E, p) = Ti[ESmip)] where Sm is the aver- 
age sequence superoperator given by Sm = J2i ^im/^^- 
Now by defining A to be the average over A^^. and by as- 
suming that SAi- = — A is small [l^l , the average 
sequence superoperator can be written as 



5,„ = Ao[Wg(A)]°™ + 0((5A), 



(3) 



where W represents the twirl over the group Q and is 
given by 



(4) 



Depending on the group Q, this map can have a simple 
structure with a small number of parameters [2^ . 

Below we explicitly show the symmetrizing effect of 
three relevant twirls. We write a map A on n qubits 
using the Liouville representation in the Pauli operator 
basis, the Pauli transfer matrix (PTM) TZ, as 



{TlAh = ^Tr[P,A(P,)] 



(5) 



where Pj is a Pauli operator for j 7^ 0, Pq is the identity 
operator /, and d = 2" is the dimension of the system 
[23| . This mapping is a matrix representation of quantum 
maps, defined over the vector space of the Pauli opera- 
tors, where composition becomes matrix multiplication. 
This allows us to write the twirl as 



We(7^A)-^^7^^,7^A7^^,. 



(6) 



The group average over the conjugation action of a group 
is well known in the representation theory of finite groups 
and, as a corollary of Shur's lemma, would be propor- 
tional to the identity if the TZ matrices were irreducible 
(i.e . there is no basis in which each TZu is block-diagonal) 
25|. In the case where the representation is reducible 
the resulting expression contains cross-terms j26j . In 
this paper we consider a bipartite system of dimension 
d = di + d2 and twirling over the following groups: C 
(the set of Clifford operators on the full space) , C (g) C 
(independent Clifford operators on two different subsys- 
tems), C (8) / (Clifford operators subsystem 1), and I (E)C 
(Clifford operators on subsystem 2). 

Full Clifford Twirl - In the case of the full Clifford 
group there are two irreducible subspaces corresponding 
to two stabilized operators: the projector onto IIq (the 
identity operator) and a projector (11) onto the remainder 
of the Pauli group. For the full Clifford group 



Wc(7^^ 



1 
an 



(7) 



which is a depolarizing channel with a = 
Tr(n7^A)/Tr(n). If we decompose the initial state 
and the measurement operator in the Pauli basis, 
p = '^jXjPild and E = '^jCjPj, (where we have 
additionally absorbed the error from the final gate into 
E), then the twirled map sequence fidelity is given by 



Fscq{m,E,p) = Aa"^ + eo, 



(8) 



Provided that ejXj is non-zero 



where A = J2j 

(we have an initial state with some polarization in the 
same direction as our measurement), the parameter a can 
be extracted by fitting the sequence fidelity to a decaying 
exponential. As shown previously the average gate 
error is r = (d — 1)(1 — a)/d. 

Cf^ Twirl - When we twirl over the group of sub- 
system Clifford operators C CS) C we obtain four distinct 
irreducible subspaces: IIq = /®/, 112 — ^®P, IIi — P®/ 
and Iii2 = P(8)P where P is the vector of Pauli operators 
for each subsystem. From Shur's lemma we obtain 



Wc»c(7^A) = 



/I \ 

a2|in2 

ai|2ni 

yO ai2lli2j 



(9) 



where ak\k' = Tr(n/;7?,A)/Tr(nfe). This is just a sum over 
tensor products of depolarizing channels, with ak.\k' rep- 
resenting the depolarizing effect on subsystem k whilst 
simultaneously twirling k' . We estimate the average 
gate error for Clifford operations in this context by 
rk\k' = {dk - 1)(1 - Oik\k')/dk [ill, where dk is the di- 
mension of subsystem k. If the errors are uncorrelated, 
A = Ai (g) A2, then the correlation coefficient is the 
product: 0,12 — Oii\2'^2\i- Consequently, any deviation, 
5a ~ ai2 — Q!i|2Q^2|1j is an indicator of correlations be- 
tween the subsystems However, some special cases 
remain undetected by this test (3^ . The twirled map 
sequence fidelity is given by 

Psoq(m, P, p) ^ Aia'{\^ + ^2a^i + ^i2a™2 + ^o, (10) 



where Ak ~ '}2jeiik ^j-'^j- '^^^ easily extract the a's 
provided we can prepare initial states (or measure oper- 
ators) with support in only one of the irreducible sub- 
spaces giving a single exponential decay. Unlike a full 
Clifford twirl, SPAM errors can have an effect, but are 
detectable by a deviation from a single exponential. 

Single Subsystem Twirl - Twirling over a single sub- 
system is more complicated. The general form is 



Wc»7(7eA) 









r 



(11) 



V 



where TIa^ ^^'^ Pauli transfer matrix of the 
map A2(p2) = Tri[A(/i ® p2)]/di and P = 
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FIG. 1. (color online), (a) General approach to characterize 
the effect of cross-talk and unwanted quantum interactions 
in a processor with multiple subsystems. The experiments 
are performed pairwise. For each pair, three experiments are 
performed: RB is performed on each subsystem while leav- 
ing the other subsystem unperturbed (b and c); then RB is 
performed simultaneously (d). 



E;^o iiPl ®h)nx {Pi ®l2)]l{dl-l). Due to the bloek 
structure of this map, repeated applications have the 
same structure, and tracing over subsystem 2 will give 
a depolarizing channel on subsystem 1 (a one parameter 
map of the form given in Eq. ([7])). However, the sequence 
fidelity is not necessarily exponential: after tracing out 
the second subsystem in both measurement and prepa- 
ration, it is given by 



Fscq{m,E,p) 



=eo + A(r")o,o, 



(12) 



where A = X^jGPigi/ ^i-^j- Fortunately, as the error goes 
to zero, the leading contribution to the 00 matrix element 
of (r™) is a'l" with ai = Tr(ni7eA)/Tr(ni). Further- 
more, if the data is well approximated by a single expo- 
nential decay, then the replacement of (r™)o,o by a™ is 
valid and the error on subsystem 1 when nothing is done 
on subsystem 2 is given by ri — {di — 1)(1 — ai)/di. We 
can obtain a similar estimate for r2 with the twirl Wc®/- 
We can then quantify the increase in error rate from 
simultaneous control. This gives us a metric of address- 
ability: the additional errors induced on subsystem k 
from controlling fc'. 



\rk - rfc|fc/ 



(13) 



This metric combines the errors from classical cross-talk 
leakage errors and quantum coupling between the sys- 
tems. 

In short, the addressability protocol is the following 
three experiments: 

Experiment 1: Implement RB on the first subsystem, 
i.e. twirling with the group C ® I (see Fig. [T}d). Fit the 
average decay of subsystem I's initial state to obtain ai 
and thus the error ri. 

Experiment 2: Perform the same experiment on sub- 
system 2 (see Fig. [Tt) yielding a2 and r2 . 

Experiment 3: Implement RB on both subsystems si- 
multaneously, i.e. twirl with C ® C, (see Fig. [TJl). Fit 



the decays of the single subsystems (e.g. ZI and IZ) 
and two-qubit correlations (e.g. ZZ) to obtain ai\2, Oi2\i 
and ai2. The addressability is then quantified from Eq. 
[Ol and potential correlations in the errors are flagged by 
examining 5a = ai2 — ai\2<^2\i- 

To demonstrate this protocol we implement it in two 
samples, a and b, with identical layouts. The sam- 
ples consist of two single-junction transmons (SJT) 
coupled by a coplanar waveguide resonator. We per- 
form single-qubit rotations by applying microwave pulses 
along individual qubit drive lines resonant with the 
qubit frequency. The pulses were shaped for high- 
fidelity control [i^]. The pulses applied to each qubit 
belong to the single-qubit Clifford group, generated by 
|X±^/2j ^±7r/2; -'^TT, ^} whcrc Rg represents a rotation 
of angle 9 around R. 

The design of all SJTs are as described in Ref. (23j . 
The SJTs on sample a{b) have transition frequencies 
a;i/27r = 4.9895 (4.7610) and lu2/2tt = 5.0554 (5.3401) 
GHz, with relaxation times of 9.7 (9.4) and = 

8.2 (9.9) ^is and coherence times of T^^ = 10.3 (7.3) 
and Tj*^ ~ 7.1 (10.2) /is. The cavity of sample a{b) res- 
onates at 7.325 (7.4269) GHz. The samples are radiation 
shielded 2l| and thermally anchored to 15 mK of a dilu- 
tion refrigerator. Sample a is the same device as used in 
Ref. (23I I where QPT performed on various single qubit 
gates revealed an average gate error of ^ 3.6%. 

Each RB experiment starts with the qubits in the 
ground state |00). In order to simulate tracing over the 
other qubit, a set of experiments is run with additional 
rotations to measure the final populations poo, poi, PiOi 
Pii . From the decay of poo + Poi and poo + Pio we simu- 
late tracing over qubit 2 and 1, respectively. The results 
of the three experiments are shown in Fig. [2] and in Ta- 
ble m We find a single exponential decay is a good fit to 
the data, with a maximum reduced of 1.63 (most fits 
have < l)i and thus the addressability RB protocol 
is valid (for more information on the fitting procedure 
see the supplementary material [1^). From these exper- 
iments the addressability error in sample a is ~ 0.5% 
whereas in sample b the error is indistinguishable from 
zero. Furthermore, from the decay of poo +P117 we can 
obtain ai2 (Table|l| and for sample a this is substantially 
different from Q;i|2«2|i which is consistent with correlated 
errors. 

We expect the qubit frequency spacing to play a large 
role in addressability. We do not have perfect classical 
addressing of the qubits; therefore, we expect the off- 
resonant drive of one qubit 's control pulse to give an AC 
Stark shift on the untargeted qubit. In addition, the 
strength of the cavity-mediated coupling is inversely pro- 
portional to the qubit-qubit detuning. In our experiment 
we expect these effects to be much stronger in sample a 
because the qubit-qubit detuning is A/2n = —66 MHz, 
whereas in sample b, A/27r = —579 MHz. This difference 
is successfully captured by our metrics. 



4 




100 200 300 400 2 64 128 192 256 320 




100 200 300 400 2 64 128 192 256 320 

Sequence length Sequence length 



FIG. 2. (color online) Experimental results of RB experiments 
on sample a (left column) sample b (right column). Circles 
(red) are individual RB fidelity decays, triangles (blue) are si- 
multaneous RB decays. Top row shows an effective projection 
into qubit f 's subspace, bottom row, into qubit 2. Fits to sin- 
gle exponential decays (solid lines) give a maximum reduced 
~ 1.63, demonstrating a reasonable fit model. The signif- 
icant increase in decay rate for sample a under simultaneous 
RB witnesses addressability errors. 



Twirl Group 


Sample a 


Sample b 




cc 


91 


0.0039 ± 0.0001 


0.0029 ± 0.0002 


r2 


/(S 


)C 


0.0067 ± 0.0002 


0.0037 ± 0.0003 




C5 


5C 


0.0086 ± 0.0003 


0.0032 ± 0.0003 




Cd 


5C 


0.0120 ± 0.0005 


0.0043 ± 0.0002 








0.0047 ± 0.0003 


0.0003 ± 0.0003 








0.0053 ± 0.0005 


0.0006 ± 0.0003 


Sa 






0.0050 ±0.0018 


0.0015 ±0.0007 



TABLE I. Summary of the extracted data from the three RB 
experiments on the two samples in terms of effective error 
rates and addressability metrics. Uncertainties are reported 
as la confidence intervals. 



Ramsey experiments conditioned on the state of the other 
qubit) and we find for sample a, mi = 0.19, TO2 = 0.32, 
/ii = -0.088, fj.2 = -0.16, 1^1 = -0.025, 1^2 = -0.048, 
and C/27r = 1.1 MHz. From these values we estimate 
(5ri|2 = 0.0034 and dr2\i = 0.007. These estimates are of 
the same order as the measured values, though it is likely 
that this simple model does not capture all addressability 
errors. We note however, that despite the higher address- 
ability error for sample a, a high two-qubit gate fidelity 
(95% — 98.5%) was measured [2J], suggesting that the 
tradeoffs between gate fidelities and addressability will 
need to be further explored. 

In conclusion, we have presented a protocol to measure 
the addressability of the gates used in a quantum infor- 
mation processing device. We apply this protocol to two 
samples of two SJT coupled by a co-planar waveguide res- 
onator. In the first sample the single qubit average gate 
error is measured to be ~ 0.5% and with operations on 
the second qubit this increases to ^ 1.0% (an address- 
ability error of ^ 0.5%) indicating that benchmarking 
a single component docs not necessarily account for its 
effect on the larger device. Furthermore, QPT on this 
device gave an even larger average single qubit gate error 
of ^ 3.6% which we attribute to the sensitivity of QPT 
to SPAM errors. The second device had an equivalent 
average gate error, but much smaller addressability error 
due to a larger detuning between the qubits. 
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J. R. Rozen, Marcus P. da Silva and Joseph Emerson. 
We acknowledge support from lARPA under contract 
W911NF-10-1-0324. All statements of fact, opinion or 
conclusions contained herein arc those of the authors and 
should not be construed as representing the official views 
or policies of the U.S. Government. 



In our system the qubits are coupled by an XX inter- 
action, with strength J, so the Hamiltonian is written in 
the dressed basis as 

H =£1 {XI + (toi2 - iyi)IX - filZX + mi2fi2XZ) 
+ £2 {IX + (msi ± V2)XI + [i2XZ - m2UixZX) 
- ijJiZI/2 - UJ2IZ/2 ± C^Z/4, 

(14) 

where £1(2) are the shaped microwave amplitude of a 
drive applied on qubit 1(2), »ni2(2i) represent spurious 
crosstalk due to stray electromagnetic coupling in the 
device circuit and package, /J.i(2) are the cross resonance 
coupling parameters (/Xi(2) = J/A for ideal qubits (28j). 
and C is an energy shift which can occur due to cou- 
plings in the higher levels of the SJTs [2^. These pa- 
rameters can be measured by a series of 7r-Rabi and 
TT-Ramsey experiments (standard single qubit Rabi and 
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PARAMETER ESTIMATION: 
ANALYSIS 



A STATISTICAL 



The experimental data is fit to a three parameter 
model using standard non-linear least squares regression. 
Wc use the nlinf it routine provided by Matlab [![. Un- 
certainties are identified from estimates of the Jacobian 
at the best-fit point with a linear assumption. To do 
this we use the Matlab routine nlparci with the confi- 
dence set to 68% (Icr) [H. Identifying goodness of fit is 
somewhat muddled because of open questions about the 
underlying distribution of sequence fidelities. There is no 
guarantee that random sequence fidelities follow a normal 
distribution. Indeed, anecdotally we observe clustering of 
high and low fidelities in some data sets. Furthermore, 
for very high sequence fidelities we see skewed distribu- 
tions because the fidelity is bounded by one. Neverthe- 
less, the residuals from fits of the average fidelities appear 
to follow a normal distribution and as such we can use a 
chi-square test for the goodness of jit. The reduced chi- 
square is given by x^/dof, where dof is the number of 
degrees of freedom and is the chi-square given by 



N 

E- 

i=l 



(1) 



Here N is the number of experimentally measured quan- 
tities Xi, fj,i is the expected value from the model under 
test and af is the expected variance in the measured 
quantities. In this case wc take N to be the number of 
truncations used, Xi is the experimentally obtained se- 
quence fidelity after averaging over the different random 
Clifford strings, and wc use the standard error in this 
estimate of the sequence fidelity as cTi . 

For sample a there are 32 different truncations and 
with the 3 fit parameters (29 degrees of freedom) we 
find for (ai, a2, ai|2, Q!2|i7 O12) the reduced chi-square are 
(1.567, 0.809, 1.369, 1.639, 1.189) respectively. The data 
for sample b have 11 different truncations which gives 8 
degrees of freedom and we find the reduced chi-square for 
these fits are (0.638, 1.609, 0.156, 0.489, 0.280) respec- 
tively. All these reduced chi-square's are close to or less 
than one so we conclude that the model is a good fit to 
the data. 



REVIEW OF GROUP THEORY 

In this section we review some of the standard results 
of the representation theory of finite groups and use this 
to decompose the PTM's of the various groups from the 
main text into their irreducible components. This decom- 
position greatly facilitates twirling (i.e. the group aver- 
age). 



Representation theory of finite groups 

Since we will extensively be using Schur's lemma later 
in this document we will give a very brief overview of 
representation theory for finite groups. A group is a set 
Q and an operation o with the following properties: 

• closure for all 51,52 G G, ,93 ~ gi ° 52 is also an 
element of Q, 

• associativity for ah 31,32,53 e G (5I 052) o 53 = 
91 ° (.92 0.93), 

• identity there exists an element I E G such that 
for all 5 G ^, / o 5 = 5 o / = 5, 

• inverse for all g E G there exists a 5"^ £ G such 



that 95 



.9 .9 



In the remainder of this document we will take the group 
operation, o, to be multiplication unless otherwise stated. 

A representation of a group is a map a from group el- 
ements to a set of complex matrices, a : G C'^^'^, such 
that the group operation is preserved under matrix multi- 
plication of the representation, i.e. cr{gi)a{g2) = (t(5i52). 
Note that multiplication can be preserved even if a is not 
a one to one mapping. For example, all groups have a 
trivial representation where cr^g) ~ 1 since cr(5i)o'(52) ~ 
(7(5152) is always satisfied by 1 x 1 = 1. A represen- 
tation that maps each group element to a unique ma- 
trix is known as a faithful representation. Additionally, 
any matrix representation of a group with non-vanishing 
determinants is equivalent under a similarity transfor- 
mation, (j{g) — > Sa{g)S^^ , to a unitary representation, 

'^{9)-'=<r{gV 0. ■ 

The final requirement before stating Schur's lemma is 
to define the notions of reducible and irreducible repre- 
sentations. A representation is reducible if there exists 



2 



a similarity transform S such that each group clement is 
block diagonal. That is, for all g £ G, 



Sa{g)S' 



ai{g) 
a2{g) 



(2) 



If there is no such S then the representation is irreducible. 
We traditionally define two representations as distinct if 
there is no similarity transformation that maps one to 
the other. Any general reducible representation can be 
written as a direct sum of its irreducible components. 
Mathematically, representations are more commonly de- 
fined by the vector space the ct's are acting on, and so 
we will generally talk of a reducible representation being 
composed of irreducible subspaces which are preserved, 
or stabilized, by multiplication by elements of the repre- 
sentation. 

There are many ways to state Schur's lemma, but a 
particularly transparent form for the follwoing is from 
Tinkham (which also contains detailed proofs) It is 
split into two components. 

Lemma 1. Schur's Lemma: The only matrices that com- 
mute with all elements of an irreducible representation of 
a group are constant matrices (i.e. scalar multiples of I). 

Lemma 2. Given two irreducible representations (Ti and 
02 o,nd a rectangular matrix M that satisfies Mai{g) ~ 
a2{g)M , for all g & G, either M = or M is invertible 
and thus oi and 02 o-fc equivalent modulo a similarity 
transformation. 

We can now discuss the implications Schur's lemma has 
to twirling in three increasingly more general corollaries. 

Corollary 3. If a is an d-dimensional irreducible repre- 
sentation of G then for any matrix M G C'*'^'' 



1 , , , 1 Tr(il/) ^ 

' ' see 



(3) 



Proof. The sum on the left hand side of Eq. commutes 
with (j{g) for g £ G- Explicitly, 

= ^ a{g')Ma {g'~^~g) , 
g'eg 

= Y,a(g')Ma {g'-')a{g), 
g'eg 

= {Y.^''^9)Ma{g)-'^ a{g), 

(4) 

where in the intermediate steps we have taken g' ~ gg. 
Due to Schur's lemma wc know that any operator that 



commutes with all elements of the irreducible represen- 
tation of a group must have the form 



(5) 



and to solve for a we take the trace of both sides yielding 
Tr(M) = ad. (6) 



□ 



A generalization of the previous corollary is 



Corollary 4. // a is a direct sum of distinct irreducible 
representations of G, CTi © (T2 ffi . . . Um then for any matrix 



i-5^a(5)Ma(g)' 

' ' geg 



j 



Tr(MP,) 



(7) 



where ¥j is the projector onto the irreducible subspace 
defined by aj . 

Proof. As in the previous case, the sum on the left hand 
side of Eq. ([7]) (which for compactness of notation we will 
label $ — '^g^g o'{g)Ma{g)~^ /\G\) commutes with a-{g) 
for all g £ G- Furthermore, 



P,a(5)$Pj =P,$a(g)Pj 
a,(5)P,$P, =P,$P,a,(5) 



(8) 



The matrix Pi$Pj is now subject to the constraints of 
Lemma [2] and thus 



(9) 



and as in Corollary |3] we can take the trace of both sides 
to obtain 



Tr(A/P,) = a, Tr(P,). 



(10) 
□ 



Finally, we can extend this corollary to the case of 
a completely general reducible representation. Up to a 
similarity transformation a general representation can be 
expressed as a direct sum of Uj copies of each of the 
irreducible representations Oj. We will explicitly write 
the subspace spanned by the /c"^ copy of representation 
Oj as being spanned by the set of vectors \vj,k.i)- With 
this convention 

Corollary 5. Given a general representation a of G then 
for any matrix M G C"^^** 



]-Y^a{g)Ma{g)-^ 



\Q\ 



see 



E E 

j k,k' = l 



(11) 



where Qj,k,k' = J2i I'^j-k.i) 



Ij.k.k' 
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Proof. This proof proceeds along the hnes of the last but 
now we use the projectors Pj.fc = J^i I'^j-k.i) {'fj.k,i\ and 
find that 

Pj,fe$Pjv,fc- = Sj,fa,,k,k'Qj,k,k', (12) 
If we multiply both sides by Qj^. j., and trace we find 

Tiiql,^,,M) = aj,k,k' Tr(Qj,^,,Q,- fc,,.). (13) 

□ 

The last ingredient from the representation theory 
of finite groups that we will require is the following. 
If we have a group G with irreducible representations 
{ci, (72, . . .}, we can form the irreducible representations 
of the group G xQ (with elements ffiffi, 5iff2, ff25i, • ■ 0^ 
the tensor products {cti (g) cti, cti (T2, cr2 (8) cti, . . .} [21 ■ 
This is of particular interest when we consider composite 
Hilbert spaces with many subsystems. 

Decomposition of some useful unitary groups 

Pauli group 

Let us consider the set of channels V with elements 
of the form Ak{p) = PkpPk and derive a general form 
for the PTM. Here, as in the text, Pk are the n qubit 

Pauh operators {/, X, Y, Z}®" with / = (^^ i)'^ = 

(? o)'^= (z o') ^^dZ= (^J \y Asacon- 
vention we will take Pq ~ /®". 

First, Pauli operators either commute or anticommutc, 
so PkPjPk ~ ±-Pj and thus 

(^aJ^. -±^Tr(P,P,) = ±'5.,- (14) 

Furthermore, by expressing mapping the four Pauli op- 
erators to a pair of bits (vw) according to {/, X, Y, Z} 
{00,01,10,11} (and therefore Pj to two binary vectors 
(v,w)) we find that for the channel above the elements 
of the PMT are 

(7^AJ,. = (-l)^*"■-''' (15) 

All the PMT's are diagonal with half the elements 1 and 
half —1, except for the identity which is mapped to all 
I's. 

This is somewhat counter-intuitive, since we have rep- 
resented the Pauli operators as a set of commuting ma- 
trices. Additionally, while the 4" matrices Pj do not 
themselves form a group under matrix multiplication the 
TlAk matrices do. The essential ingredient is that the set 
of channels V forms a group under composition because 
the conjugation action of a unitary channel makes them 



insensitive to a global phase (i.e. it is irrelevant whether 
X times Y is equal to ±iZ for channels). 

We can now explore the structure of the TZa^. repre- 
sentation of P in a little more detail. Group multipli- 
cation TZa^ = TZAiTi-Aj is described by (v'*^' , wC^'^) = 
(v'^^'' v^*'', w*^^^ ® w^*)), and is simply a representation 
of Zj®^". In fact, the TZa^ matrices are what is known 
as the regular representation of l®"^"^ which is defined as 
the representation over the vector space of the group el- 
ements themselves. For 1®"^"^ the regular representation 
is the direct sum of all 2^" distinct irreducible represen- 
tations. 

Clifford group 

The Clifford group C is the normalizer of the Pauli 
group in the space of unitary channels, that is, conjuga- 
tion by an element of the Clifford group at most permutes 
the elements of the Pauli group or CjP^C] = ±Pk for all 
Pauli operators Pi. Therefore, the PTM's describing ele- 
ments of the Clifford group have a 1 or —1 in each row and 
column with the remaining elements zero. Furthermore, 
conjugation preserves commutator relations, or alterna- 
tively quantum channels are completely positive, which 
places additional restrictions on the set of physical TZ 
matrices. 

We will now show that the PTM representation of the 
Clifford group is of the form I © cr, i.e. a direct sum 
of the trivial representation of the Clifford group and 
an irreducible representation we will label a. The trivial 
irreducible subspace corresponds to the identity (Pq = /) 
which is left invariant by all Clifford channels. To show 
that the remaining subspace is irreducible we note the the 
P is a subgroup of C. The implication from the previous 
section is that the only matrices that commute with all 
elements of the representation are diagonal. The only 
restriction on how the Clifford group permutes the Pauli 
group is that it preserves commutator relations so for any 
two Pauli matrices Pi and Pj ^ 0)there exists an 
element of the Clifford group, C, such that CP.C'' = Pj. 
Therefore no subspace of the Pauli's is stabilized except 
for the entire group (minus the identity). 

For a composition of Clifford operators acting inden- 
pendently on n-subsystems, C®", the irreducible repre- 
sentations are of the form (I -|- cr)". As an example, 
the group C (E)C has the following irreducible representa- 
tions: II = {//}, 1(7 = {IX,IY,IZ}, al = {XI,YI,ZI} 
and aa = {XX, XY, X Z,YX,YY,Y Z, ZX, ZY, Z Z}. 
Addionally we could consider the Cliford group act- 
ing on a subsystem of a large space. In this case 
we are still representing the single subsystem Clif- 
ford group, but with degeneracy. For example, 
I (E) C has four copies of the irreducible representa- 
tion I ({//}, and {ZI}) and four copies 
of cr {{IX,IY,IZ},{XX,XY,XZ},{YX,YY,YZ} and 
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{ZX, ZY,ZZ}). 



Unitary group 



The structure of the full unitary group is essentially 
identical to that of the Clifford group, that is the PTM 
representation of a general unitary channel has the form 
I © (T, where the identity is stabilized and the remaining 
subspace is irreducible. This should come as no surprise 
since the PTM mapping is quadratic in the unitary ma- 
trices and the Clifford group is a 2-design of the unitary 
group. In fact, requiring that all but the identity forms 
an irreducible subspace is equivalent to the statement 
that the Clifford group is a 2-design Q. 



Twirling 

Twirling a channel A over a group of unitary channels 
Q = {Ui,U2, . ■ .} is described by 



which in the TZ representation is 



1^1 u% 



(16) 



(17) 



With a detailed understanding of the representation 
structure of the Pauli transfer matrices the results of 
twirling are an almost trivial consequence of Corollaries 
[31141 and [51 



Pauli twirling 

The PTM representation of the Pauli group is a direct 
sum of 4" 1-dimensional distinct irreducible representa- 
tions (i.e. the Pauli operators themselves are each sta- 
bilized by each element of the Pauli group in this repre- 
sentation since the group is Abelian). Therefore twirling 
over the Pauli group is a direct consequence of Corollary 
2 and is described 



(18) 



Only the diagonal elements of TZ remain. This is exactly 
the space of Pauli channels, that is channels of the form 



Clifford twirling 

The Clifford group has a representation that is the sum 
of two irreducible blocks, the identity and everything else. 



For any trace-preserving map the element of TZ mapping 
the identity to itself is unity and therefore we can use 
Corollary 2 again to find 



1 if i = 

(Tr(7^A) - l)/(d^ - 1) otherwise. 

(19) 

This is the set of depolarizing channels, Adop(p) = ctp + 
(1 — a)I where the survival probability a = 4^^. ^ ^ ■ 



Subsytem Clifford twirling 

Here we will focus on a two qubit system and look at 
twirling over the groups C (g) C, C ® I and I ® C. The 
simplest is C ® C which we previously determined has a 
PTM that is a direct sum of four disntinct irreducible 
representations: II,Icr, cri and aa. If we define Pn, Pfc, 
PctI and P^o- as the projectors onto the respective irre- 
ducible subspace, than 



1 



Wc®c(i?A)^ 



Twirling over C C 



Tr(i?(A)Pi,)/3 
Tr(i?(A)p^j)/3 

Tr(i?(A)p,,)/9 



if 
if 
if 



P^ e {//} 
e {IXJYJZ} 

P., e {XI, Yi, zi} 

otherwise. 

(20) 

) / is quite a bit more complicated due 
the multiplicity of the distinct irreducible representation 
in the PTM. In keeping with the conventions of the first 
section of this supplement let's label the irreducible rep- 
resentations of the PTM of C (x) / as 

aij ={//}, = {IX},ai,Y = {IY},ai,z = {IZ}, 
a, J ={XI, YI, ZI}, a,,x = {XX, YX, ZX}, 
G,^Y ={XY, YY, ZY}, G,^z = {XZ, YZ, ZZ}. 

(21) 

While we could derive twirling from Corollary 5, there is 
a more intuitive way to calculate this expression which 
we will show here. 

First we can express the two qubit TZ matrix in terms of 
the single qubit Pauli operators by doubling the indices, 

(7^A)y■,fc^ = jTr[P, ® PjK{Pk ® Pi)]. (22) 
and therefore the twirl is of the form 



ij,kl 



1 

4^ 



E 



Tl■[UmP^Ul®PJA{UmPkUlm)]. 

(23) 



From this expression we can show that if i 7^ fc the ma- 
trix element is always zero by a symmetry arguement. 
For any two Pauli operators Pi ^ Pk (with Pk ^ Pq for 
this example) there is always an element, C, of the Clif- 
ford group such that CP^C^ = P^ and CPfcC^ = -Pfc. 
This is due to the fact that the only constraint on 
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the Clifford group is that it preserves commutation re- 
lations. Substituting Um UmC in Eq. ((23)) yields 



Wc®/(7^A)ij,fe/ = -Wc®/(7^A)iJ^fc/ = for i ^ j. For 



the non-zero components we note that the Clifford group 
preserves the ideneity and unfirmly redistributes X, Y 
and Z so that the full form of the twirl is 



The interpretation of this quantity is that for each con- 
figuration of the untwirlcd qubit we have a single qubit 
twirl. 



Twirling over the entire unitary group yields exactly 
the same result as the Clifford group. 
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